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ABSTRACT 


This report is intended as a rigorous tutorial on synthetic aperture radar 
(SAR) with emphasis on digital data collection and processing. Background 
information on waveform frequency and phase notation, mixing, I, C conversion 
sampling and cross-correlation operations is included for clarity. The fate of 
a SAR signal from transmission to processed image is traced in detail, using 
the model of a single bright point target against a dark background. Finally, 
some of the principal problems connected with SAR processing are discussed. 
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SYNTHETIC APPiRTURE RADAR AND DIGITAl, PROCESSING: 

AN INTRODUCTION 


I. Introduction and Motivation 

The purpose of this report is to provide a basic yet rigorous 
introduction to synthetic aperture radar (SAR), with an emphasis on digital 
data collection and processing. The method will be to start with simple models, 
and then to refine these models so they more closely approximate the real 
world. 

The only knowledge assumed on the part of the reader is a familiarity 
with elementary algebra, trigonometry, and complex arithmetic, plus a faint 
glimmer of differential calculus. All other concepts will be developed in this 
report. 

Completeness is not attempted here. Instead, we trace in detail what 
happens from beginning to end when a SAR interrogates a point target on the 
ground. The principle of linear superposition can then be used to infer the 
result when th.^ terrain contains distributed targets, i. e. . multiple point 
targets. In the last section, some problems connected with SAR are sketched 
for the interested reader. Some knowledge of probability is assumed. 

We begin in this section by describing a simple kind of imaging radar; a 
£e^ aperture radar with no rang ' compression. While this radar is not a SAR, 
an understanding of how it generates imagery will prove valuable in the under- 
standing of a SAR. 

Our primitive radar may be thought of as a flashlight pointing downward and 
outward from the side of a moving airplane. Figure I shows a view of the beam as 
seen from directly behind the aircraft. The angle is called the "range" ele- 
vation beamwidth and is related to the grenmd swath width. Of course, the beam 
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Fig. 1. Radar Beam as Viewed from Behind the Airplane 

is composed not of light, but of lower frequency radio waves, and the 
"flashlight" is really a radar antenna. Figure 2 shows the beam as viewed from 
directly above the airplane. Notice that we have drawn the beam considerably 
narrower in Fig. 2. The small angle is the azimuth beamwidth . With the 
proper design of a radar antenna, it is indeed possible to create a beam with 
distinct range and azimuth beamwidths. 

Now let us suppose the radar is turned on for a very short time, and then 
off again. A brief burst of radio energy is emitted, which we may picture as 
the shaded band in Fig. 3. (Note incidentally in Fig. 3 that ground range G is 
related to slant range R by 
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R ^ 












where is die angle of elevation of the target from the nadir. After reflection 
from the ground, the same shaded band is shown returning to the aircraft in 
Fig. 4. Note diat only the two designated areas (X) on the ground have enough 
reflectivity to return an appreciable echo, and that the two echoes will arrive 
back at die aircraft at different times, because diey must travel different 
distances. Suppose we have a receiver which samples die returning echo 
as a function of time and records its instantaneous intensity. Then the 
sampled intensity will appear as in Fig. 5, or as in Fig. 6. Thus we see 
that by sampling the echo received from a single burst of radar energy, we 
have ’’painted” a one-dimensional strip of imagery perpendicular to the air- 
craft flight path. (This ’’strip” will henceforth be called a ” range line” in this 
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report, since it is oriented along the rang e direction, i. e. . perpendicular to 
the flight path. The term range line will also be used later for any strip of 
unprocessed or processed radar data which is oriented along the range axis. 

The opposite notion is an azimuth line . which is a slice of data oriented in a 
direction pa rallel to the flight path. ) The paintbrush mo%'es from alniost beneath 
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Fig. 5. Sampled Plot of Received Echo Intensity Versus Time 
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Fig. 6. Alternate Echo Display 


the aircraft to far out to the side; the "movement of the paintbrush arise, frona 
the different arrival times of the echo from different portions of the swath. 

Note in Fig. 6 that it the two reflectors (X) were closer together than the 
width w of the pulse in Figs. 3 and 4. the two targets would blur together. 

Thus it is imperative to transmit a very short pulse. 

Now we have a strip, as in Fig. 6. but how do we produce a two- 
dimensionai image? The answer it easy: just keep painting strip, next to 
each other. We do this by turning on the radar again, when the aircraft has 
moved over, to produce the neighboring strip. Note that the "width- of our 
strip in azimuth (Fig. 6) is equal to the azimuth beamwidth along the ground. 

We cannot "see" any detail within the bright spots of Fig. 6 because the 
sampler can only record one intensity value at each Instant of time. Thus the 
two X's in Fig. 2 cannot be separated. If we now continue, turning the radar 
on and off, we will "paint" a long succession of strip, next to each other. The 
azimuth width of each strip is the ground azimuth beamwidth; and in the imagery 
we can only distinguish objects in one strip (i. e. , one pulse) if they are furt 
apart in range than the range pul.ewldth w (actually 1/2 the pulsewldth - «n you 
see why? ). Thus our image resemble. Fig. 7, if we allow four "strips. " 
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There are problems connected with this imaging approach. The first is 


that the pulse widtli in range (w in Fig. 3) must be kept short, as stated above, 

Tliis requirement limits the energy transmitted per pulse and makes it difficult 

to detect weak, distant targets. The second problem is that, as already stated, 

we cannot distinguish tlie two targets (M) in Fig. 2. Although we can alleviate 

this problem somewhat by making |3 , the azimuth bcamwidth. very narrow, we 

a 

can only succeed partially, since very large antennas are required to produce 
narrow beams. Synthetic aperture radar (SAR) is a means of solving tlie 
sec ond problem. In short, a SAR is simply a radar which uses the information 
from an entire sequenc e of pulses along the a/.imuth (and not Just one) to produce 
single painted strip of imagery as in Hg. h. Thus, for example, all 4 pulses 
implied by I'ig. 7 miglU be used in a SAR to provide infortnation al)o\it the first 




output strip in the figure. The result is refined resolution or ability to 


separate nearby targets along the azimuth. 

The first problem (i, e. , the problem of weak signal returns due to a 
short duration pulse) is solved by simply transmitting a longer pulse, often of 
the type known as the "linear FM chirp. " The reader may well ask at this 
point, "Then how do we obtain the required range resolution? " The answer is 
that upon receiving the echo pulse we " compress " it to a shorter " pulse" by the 
technique of cross-correlation, to be discussed in section III-B. However, 
before discussing cross correlation we first require some background on 
electromagnetic waveforms. The reader is urged to follow the detour method- 
ically; the results will be well worth the effort. 


II. 


KloctromaKiictit U a\ cl'orm Notation and Discussion 


A. Phase and Frc(|Uoiu y 

Radio waves, like light waves, travel at a speed given by 

c ' 2. 9‘^77o X 10^ mcters/scc, (1) 

and are typically expressed as sinusoids: 

s(t) = A cos (u^t) (2) 


Here 

t = time, a real number 
A = amplitude factor, greater than 0 
u) = radian frequency, a real number 

The right hand term in parentheses, (wt). is the argunu'nt of the sinusoid and 
is called the phase. If we put, for the sinusoid of Pq. (2), 

(^(t)=^t ( >) 

then we have 


cV(t)=f 


(4) 


Thus tho rciclian Iroquonuy is the tinu' durivati\ o of the phase. In other words. 
frec|ueney is tlie rate of etuu^ce of the ph.ise. ^ 

A more j;oneral wa\eform is given by 


s(t) - A(t) eos (6(t)) 


(5) 


0 


whore both A and 6 arc real functions of t. Again, <1) is called the phase 
fum tion and the radian frequency. The (positive) function A(t) is the 

amplitude function. If we change variable by forming a new function f(t) 

4>'(t) = 2TTf(t) (6) 

then f is the frequency in Hertz , or cycles per second. Thus in the special case 
of Eq. (2), we have 

4>(t) = cot 

«t>'(t) = cj, a constant function of time. 


We have then 


f(t) = co/2it = f, a constant 


(7) 


so that 


s(t) = A cos (27Tit) (8) 

From (8) wo see that s(t) has period r, where 

^=l/f (9) 

Example 

Tlie following waveform is known as a linear FM waveform. ("Linear FM" 
stands for linear frequency modulation. ) This form of signal, also called a 
"cliirp, " is frequently used in SAR's. Define, for -T/2 ^ t T/2 and real a. 


l.(t) 



(10) 


10 



Note that at t = 0 we have 


f(t) = 4>'(t) = 0 


wliile at t = we have 


4>'(t) = ± 


aT 

2 


(14) 


(15) 


Note also that we may differentiate the frequency in Eq, (12) and (13) to obtain 
the frequency rate ; 

<b"(t) = a ( radians/ sec^) (16) 

f (t) = f'(t) = ^ (Hertz/sec) (17) 


Eq. (16) gives the slope of the line in Fig. 8. 

B. Interference and Antenna Beamy, Idth 

Suppose we have two waveforms of unit amplitude and unit frequency in 
Hertz: 


Sj (t) = cos (2TTt) (18) 

s, (t) = cos (ztt (t - •l-lj (19) 

See Fig. 9. We see that S| and s, differ only in phase, but this is enough to 
cause total cancellation. Thus, if we put 


we base 


s(t) = Sj (t) 1 - sjt) 


s(t) 0 


( 20 ) 


( 21 ) 


) 


Fig. 9. Interference 

We note from Eqs. (18) and (19) that Eq. (21) holds when the offset 
bctNveen the signals is one-half of a period. If we regard a wave as a physical 
structure, then we may define a wavelength as the physical distance corre- 
sponding to "the smallest period". We usually denote this quantity by and 
observe that tlie "physical wavelength" \ and the "time wavelength" t (also 
called tlie "period") are related by 


1 


where c, the velocity of light, is given by (1). Eq. (22) is only a version of 
the familiar relation 

distance = (rate) x (time) 

and if we recall that period is inversely related to frequency by 

T=i (9) 

f 

we obtain from Eq. (22) 


In words, Eq. (23) says the wavelength times the number of waves passing 
per second is equal to the distance light travels in a second. 

If in Fig. 9 we regard Sj and s, as signals possessing the same 
arbitrary wavelength then we see that in general tv-o signals will interfere 
destructively ("cancel") when they are offset by X/2. 

Antenna beamwidth is related to the above considerations. In f'ig. 10 we 
show a point target returning a radar echo (of wavelength \) to a radar antenna 
of length L. The target is offset from the perpendicular axis ("boresight") of 
the antenna by 0 radians, and is sufficiently far away that it produces a nearly 
straight wavefront of parallel waves at the antenna. In the figure we suppose 
that the angle 0 is such that the wavefront is \ units further from the right 
edge of the antenna than from the left, (i. e. . the offset between waves and 

is a wavelength). The offset between waves and is \/2. and 

they cancel out. Looking, for example, at Q and @ we see that every 
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/ / BORESIGHT / 

WAVES FROM DISTANT POINT 
TARGET 

Fig, 10. Destructive Interference in a Radar Antenna 

wave in the left half of the antenna has a corresponding cancelling wave in the 
right half, and thus we expect the net signal to be zero. This is indeed the 
case, and if we plot the strength of received signal from a point target as a 
function of offset angle, we obtain a curve resembling that in Fig. 11, where 
G satisfies (see Fig. 10) 

sin 6 ' (2 

Since generally in SAR antennas, \ is much smaller than L. 

\ cc 1, 


15 



Fig. 11 . Antenna Pattern 


.U sin e ^ G and hence Eq. ( 24 ) beconr.es 
e have that 6 is small, sin 6 - « 


, f lO for spaceborne SAR systems, 

CeneraUy th. be.n,widtH e U on the order o. 1 .P 

so ( 25 ) is a very good approximation, ) 

Thus the nuU-to-nuii beamwidth o£ the antenna ts 

20 =^ 

- 

maximum ,tven as 2(0,2, = 0. henee the 

"half-poNver’ or cit^ 


Ath 


1— *-a 


2 


and Nve see that the longer the physical antenna, the smaller the beaimvidth. 

Thus in Section I, where we desired a narrow azimuth beamwidth we see 
that a SAR antenna with long azimuth dimension was required. 

How long? Well, if we wish to separate Uvo targets 10 meters apart in 
azimuth, then our "flashlight beam" must be narrower than 10 meters on the 
ground along the azimuth. (See Fig. 2. ) A typical radar wavelength \ might be 
0. 235 meters, and the distance between airplane and ground target might be 
15 km. Using the relation arc length = (radians) X (radius), 



we obtain 


lOm = 0 X (15 km) 


so that 


0 = 0. 00067 radians 


and from Eq. (2b) 


L 



L = 352. 5m 
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Most airplanes are no. equipped to carry an antenna this longl A. wo have 

already stated, SAR provides a solution to the prohlent; however, our hac - 

. In the meantime note that 

around development must firs, be completed. 

Eq. (26) implies that asimuth resolution 6^ is given by 

R\ 


^A= L 


{ 26 ') 


where R is the distance (slant range) between aircraft and target. (See Fxg. 3. ) 
The resolution of a SAR will be considerably better. 

^ers alter the frequency of a sinusoidal signal by mlxln,. which con- 
sists merely of multiplying the given signal by another signal and 
In this way a signal having the sum or difference frequency iS pro uce . 

our application we require the difference frequency. 

AS an illustrative example, consider two signals having frequencies o 
HS and 120 Ha. If these signals are multiplied together, we obtain a new 
a,„ai consisting of two components, or ..sidebands. The first component has^ 
frequency (120-100) or 20 Ha. The second component has frequency 
ac 220 Ha. A lowpass filter will Isolate the 20 Ha component. 

TO understand how mixing worhs, we consider the following well Rnown 

trigonometric identities* 


cos(a + b) = cos a cos b - sin a sin b 

(27) 

cos(a - b) = cos a cos b + sin a sin b 

(28) 

sin(a + ta) = sin a cos b + cos a sin b 

(29) 

sin(a - b) = sin a cos b - cos a sin b 

(30) 
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Adding Eqs, (27) and (28) and dividing by 2 gives 


cos a cos b = j |cos(a + b) + cos(a - b)) 


( 31 ) 


Subtracting Eq. (27) from Eq, (28) and dividing by 2 gives 


sin a sin t> = ^ (cos(a - b) - cos(a + b)) 


(32) 


Adding Eqs. (29) and (30) and dividing by 2 gives 


sin a cos 1> = ■!■ (sin(a + b) + sin(a - b)) 


(33) 


From Eqs. (31) to (33) we see that multiplication of sinusoids translates to the 
summing and differencing of their arguments. 

To see how Eqs. ( 31) to (3 3) niay be used, suppose a radar is receiving 
an echo of the form 


e( t) = cos t) (34) 

where 

w, = constant freqviency 
t = time 

L-et us split the return signal into two cl'.annels, as in Fig. 12. W'e 
multiply the top cliannel signal l)y cos(., ,t) '"'d the bottom i hannel signal by 
sin(w^,t), where this time is a constant frequency of our own choosing. 



1/2 [sin (u)^ + 

+ sin - W 2 H] 

sin (^ 2 ^) 

Fig, 12. The First Stage of Mixing 

As is seen from Eq. (31) and Eq. (33), each channel output is a sum of 
sinusoids of high and low frequencies: 

cos |^(u)^ + )tj + ^ cos 
for the top channel, and 

j sin |(wj +‘^ 2 ^^] T ‘ 

for tlie bottom channel. 

L.et us suppose that - w^) is the low (positive) frequency. Then if w'e 
low pass filter each channel we obtain the output in Fig. 13, Normally, the 
term mixing refers to the use of the top channel, but if we use both channels as 
in Figs. 12 and 1 >, tlie net result is as sb.own in Fig. 14, and consists of two 
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LOWPASS 

FILTER 


1/2 co» (u), + Wj)t *■ 1/2 cot (w, - ujj)t 
TOP CHANNEL 


1/2 coi [(w, - Wj)t] 


BOTTOM CHANNEL 

1/2 tin (u)^ + tjj)t + 1/2 tin (u>^ - tjj)t 


1/2 tin[(w, - Wj)t] 



LOWPASS 

FILTER 


Fig. 13. The Second Stage of Mixing - Lowpass Filtering 


cot (Wjt) 



Fig. 14. Offset to I, Q Conversion, Difference Frequency 
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sop.n-atv outpviU. k.unv.. as in-phaso (1) and quadraturo (Q). Noto ll.al in 
Fip. 14 tl\i' scaling (ainpliUnk') constants have been removed. 

The ordered pair ot real numbers 


jeos |^(wo, - 


I 

is usually abbreviated by mathematicians as 

j [(«.', 

I* 

whore 

j = N -T 

o = the oxpoJiontial liuu tion 

Tluis the net result of mixing with a sine and cosine channel is shown in I ig, 15 

as a complex output signal witli a new frequency. 

Note that it wouhl be just as easy to use a highpass filter instead of a 
lowpass filter; in this case we would obtain tbe sum frequency instead of the 
difference frequency a.ul Figs. 14 and 15 would be replaced by Figs. In and 17. 

We are now ready to see how mixing is use.l to preserve coherency (i. e. . 
a/imuth phase information) in a SAU system. For tlie transmitte.l pulse we use 
a linear FM chirp modulating a high frequency carrier wave. 


I 


In general, for real 0. we put 


= eos 0 I i sin t) 


whit h is formally equiv.dent to a cartesian pair 


ih t • ,J 

e- - 'cos 0. sin t'^ 


nrovitle.l we operate on the components properly in subsequent 
uUUo,.., Nou. tlMl tl»- ( . ) .i.n in . 0 . ,, »in 0 i. pni-oly n torn..-.l,ty; .U 
\nd imminarv numbers ' do not Ideml. " We In.ve, tor example ' ' 

" .,r*inyrn 4,,! ...nl 1. ..qnivnl..,,i to ( ^ 41. Tl..- ..quivaUMu- .■ va..,l prnvi.U-.l 

w e (iefine multiplit ation by (.i.l>Ht . d) - (at - bd. ail I In ). 
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0 

Fig, 18. Linear FM on a Carrier - The Frequency Plot 

By way of example, it may be helpful to know that if we define tlie carrier 
frequency f^^ in Hz by 

then a typical value of f^^ for a real radar might be 1275 MHz. A typical value 
^max 1284 MHz; might be 12ut) MHz, The signal bandwidth 

would then be f - f = 1284 - I2u6 = 18 h'Hz. 
iiiax rnin 

To continue, we require that even when Lq. (38) is not being transmitted, 
a stable local oscillator (STALO) continuously generates a carrier waveform of 

freqviency = 2n f^^. This waveform is generated in two channels, producing 
signals Cj (t) and CqII). where 



Cj(t) = C 08 (w^,t) 

(3ba) 


C^lt) = sin(^^,t) 

(3bb) 

for all time t. 





In a SAR transmitter, tlu' STAl.O must lu' extremely st.ible beeause the 
carrier wave is mixed with the returning radar echo. This ntixing provides the 
critic. »l ".\zimvitl» ph.ise factor" <l' ^ in tlie tollowiim way, (Tlie use ot the azimuth 
pi\ase factor is given i«i section IV. ) 

Suppose the t< *rain is dark except for one bright point target PT on the 
ground at a liistance R meters from the radar transmitter. (This distance is 
called slant range, see Pig. 3. ) It takes tl\e signal in Ecj. (38) a certain titiie i3t 
to traverse the round-trip distance 2R to PT and back. (We make tlie sin'tplifying 
assumption tliat the aircraft does not move between transmission and subsequent 
echo reception: thus die radar signal travels R units in both directions. ) Since, 
at the speed of light c, this time of traversal is given by the relation 

c(C.t) = 2R 


WO have 



o 


(40) 


riuis the leading edge of the pulse, wliicli returns M seconds 

r r iw 

mission, is roooivod at timo f M = - “ f • Similarly 

r ' R 

of the pvilse returns at t — ♦ ^ . after wl\ich time tlie radar 
riuis in tlie interval 


.ifter its trans- 
the trailing edge 
receives no eciio. 


1 ' 






(41) 


tl\e radar receives the echo c(t), whicli is the transmitted sign.d delayed by 


elt) 


As (t 



A 





(42) 


M soioiuis: 


(The letter A denotes a positive amplitude factor proportional to the target's 
reflectivity, or "cross section,") (Note again that the correctness of (42) 
depends upon a simplifying assumption: we assume that R does not vary with 
t during reception of a single pulse echo. Without this simplifying assumption, 
the range signal mixing would be complicated by a "doppler offset frequency" 
which is normally neglected in SAR processing. ) The returning echo is divided 
into two channels (I and Q) for separate mixing with the carrier components 
Cj(t) and C^(t), The results are indicated in Fig, 19. 

As before, we represent the pair of outpvits I and O from Fig, 19 as a 
complex osponential e^(t): (the subscript r denotes range signal, since the 


e (t> = A oot ^ w I (t - j 



1 ig, I", Result of Mixing tlio Returning Radar I a lio W itii 
tlu' ST.M.O Signals Cj and 



independent variable t is "oriented" along the "range" direction. Note that e^ 
is the return from a single pulse). 


e^(t) s A exp 




Using the relatione^''’*’ = e^ e^ we may rewrite Eq. (43) as 


e^(t) s A exp 


I . I l a/, 2 r\ 


where 


n- |.ral 


is the azimuth phase factor and 


<t>. 


2.R 


( 43 ) 


(44) 


(45) 


(46) 


is its phase. We remind the reader that Eq. (44) represents the demodulated 
echo only on the interval 


T 



+ 


2R 


c 


Outside this interval e^(t) is identically zero. 

We make first a few observations about Eqs. (44) to (4u). First, the 
azimuth phase factor is independent of time t of reception. As long as the 
receiver is receiving signal from the point target PT. the phase factor Eq, (45) 
will appear as a constant term witli phase given by Eq. (4n), Next, the azimuth 




phase in L’q. (4o) may be rewritten in terms ol wavelength \. To see this, 
note that 


4> 


a 


R 2tt£ R 

o_ = . 2 X — 

e c 


and so 


4> 

a 


4nR 

\ 


(47) 


where we have used the relation = c (cf. (23)). Finally, the factor 
known *\is tho rani;o phase taetor and wc denote its phase by : 


i> - 
r 


Note tliat is ol tlie form 


(48) 


(49) 


j. - 

O - — V 


(50) 


w here 


*ind 


y = t - ^ = t - ^ 


(51) 


)' = 


2R 


(52 


IS a function only of target separation R from the radar. 




D, Samplint’ 

Lip to this point \vc have assumed tliat all signals have been analog. 
Actually, sampling can take place at a variety of points within the system; 
we choose to introduce it right here for simplicity. Thus we assume that the 
analog signal Eq. (44) is sampled as it is produced. 

The subject of sampling is well treated in digital signal processing texts. 
The main issue is sampling frequently enough so that all the information in the 
analog signal is retained. Such ideal sampling is possible only foi certain 
types of signals. One svich type is a signal having a contimious Fourier 
transform of finite extent along the freqviency axis: the interested reader may 
find the corresponding "Nyquist sampling theorem" in [11 pp. 26-30. Another 
such signal is one whose Fourier series has only a finite number of terms; a 
version of the Nyquist theorem for this second type of signal is found in (2). 
pp. 25-28. 

The gist of these theorem.s is that w e can keep most of the analog infor- 
mation of Eq. (44) if we sample it at sampling frequency f . where 

s 



Referring to Fig. 18 we see that f^ is the vertical extent of the frequency 
graph, expressed in Hertz. Thus we are saying we must sample the complex 
echo Cj,(t) at a frequency equal to its "bandwioth. " For further discussion of 
sampling theory and alternate definitions of bandwidth, the reader is referred 

to in. ui. Ml. 


The sampling period At is defined as the time interval between adjacent 
samples, it is the reciprocal of hence 




( 54 ) 
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(Normally we sample at a slightly higher frequency than the bandvvidth, and At 
is correspondingly smaller. ) 

We now reexamine the situation thus far. We have followed the radar sig- 
nal from the transmitter to a point target and back again, and seen how the sig- 
nal is mixed down (demodulated) to baseband. ("Baseband" implies the signal 
is centered about zero frequency; the carrier frequency was removed by the 
mixing of (44)), Next the signal was digitized (sampled) at intervals of 
At seconds. All of this manipulation is concerned with the return from a single 

transmission or puls e. If we knew the terrain had only one point target at 

T 2 R 

range R, we could begin sampling e^(t) at --j + — and stop sampling at 
T 2R 

+ ~ . The sampling window would then be T seconds wide, the width of 

the range chirp. However, targets at a variety of slant ranges will be return- 
ing radar echoes; thus we must increase the size of our sampling window to 
accommodate returning pulse echoes from across the entire ground swath illu- 
n^inated by the antenna in the cross-track direction. Thus we may assume we 

start sampling at tinne T^ . We take a new sample every *^t seconds until time 

T ^ R 

T, . at which time we tvirn off the sampler. We assume T v --r + — and 
“ I 2 c 

T 2R 

*2 ^ T, , The time interval over which we sample is thus AT = T, " * 

the number of samples is N where 



(55) 


Time T^ is the time ot reception ol the leading edge of the echo from the inmost 
point on the swath, an<l T, is the time ot reception of tlie tail of the echo from 
tlu’ outermost point on tlie sw atli. 
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If PT at slant range R is the only reflecting target, then our array of N 


samples looks like Fig. 20. 

The first 0 in Fig. 20 corresponds to t = . The last 0 corresponds to 

t = T, . There are N entries; hence we may denote them as a vector 
V = jv(l), v(2), .... v(N)j. Let us suppose the kth element of v is the beginning 

of the echo from PT: v(k) ~ ; we also choose 2 so that 

( T 2R\ 

~ I corresponds to the end of the echo return from PT. Then 
we have v(n) = 0 for 1 < n < k and v(n) = 0 for k + f < n < N. provided PT is the 


only reflecting target. Note that v(k), v(k+l) v(k+f) all have signal data 

from PT. How can we "compress" v(k) v(k+0 into one or just a few 

samples so as to obtain fine range resolution, as promised earlier This com- 
pression, it turns out, is the function of the cross -correlation process, and is 
discussed in the next section. 


(k - 1) ELEMENTS 
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N ELEMENTS 





NOTE: IN THF ABOVE VECTOR, ^ 


Fig. 20. Array of Samples from a Sit.gle I’ulse Echo from a 
Single I’oinLTarget 


Ill, Vectors. Matrices. Cross-Correlation and the Ranee Processing 

In this section we develop the ide&s of vectors, matrices and cross- 
correlation. While the reader with sufficient background may wish to skip over 
the explanation on vectors and matrices, the section on cross-correlation should 
be read carefully. Facts about cross-correlation that are especially relevant 
to SAR processing will be emphasized. 

A. Vectors and Matrices 

A vector , as described in the previous section, is simply a list of numbers 
in a particular order. Thus (1, 0, -3) is a vector having three components. It 
is distinct from the vector (-3,0, 1). The vector v from the previous section 
has N complex components: 

v = (v{l), v(2) v(N)) (56) 

It is customary to denote vectors as in Eq. (5b); the parentheses enclose the 
whole vector and separate parentheses are used to count the elements of the 
vector. Thus v(7) is the 7th element of 

Vectors may be added if they have the same length. Thus if 

u= (u(l) u(N)) and V = (v(l) v(N)j then we define the vector u + v 

by 

(u + v)(k) = u(k) + v(k) for k = 1, 2, .... N 

Jvist as naturally we may multiply a vector by a scalar (a scalar is just a 

number). If u is a scalar and u = (u(l) u(N)j is a v'ector then ttie 

vector axi is defined by 

(cm)(k) - a X u(k) for k = 1, 2 N . 
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A matrix or array is just a list of vectors of the same length. Typically 
the elements of a matrix are identified by affixing row and column subscripts; 
thus m^ might be used to denote the element in the ith row and jth column of a 
matrix M. For example we might define 



so that m^j = -1 and m^^ “ 3. 

The reader will note that if the digitized echo return from one radar 
pulse forms a vector v of length N, then the collection of returns from K pulses 
forms an array of N X K elements. We will (arbitrarily) assume that each 
pulse is represented as a (vertical) column in the matrix, so that the row index 
is related to the time of reception of a sample from a single pulse, while the 
column index identifies the particular pulse in question. 

B. Cross-Correlation and the Range Processing 
Let 


u = (u(l), u(2) u( m)) 


and 


V = (v(l), v(2) v(m)j 

be two vectors of the same length m. We define their cross-correlation vec'.ux 
C^^, also of length m, by 

m 

CuvC<) = ^ u(k’) v(k + ^) for k = 1 , 2, . , . , m (S7) 

1 = 1 
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where 


_ (n-1) mod m for n ?! m 

n means < 

m for n = m 

so that 

(k+^-1) mod m for k+^-1 i m 

kHC =<^ 

m for k+ f- 1 = m 
For example, if 


u = (1, -1, 0, 0, 0) 
V = (0, 0, 1, -1, 0) 


then 


=1X0 + (-1) XO + OXl + 0X (-1) + 0X0 = 0 
= lX0 + (-l)Xl + OX(-1) + OXO + OXO = -1 


^uv< = 0x0 + OxO+lXl+(-l)x(.l) + 0X0 = 2 


(58) 


= lX(l) + (-l)XO + OXO + OXO + OXl = -l 


= 1 x0 + (-l)x0 + 0X0 + 0Xl + 0X(-1) = 0 


In other words, to cross-correlate u and v. we successively slide v 'around 
itself" and then multiply by u and add. as in I'iy, 21. 

Note that has a maximum at 3), indicating that the third sliding 
position of v gives the "best match" with u. 


u 


( 1 ,- 1 , 0 , 0 , 0 ) 


U 4 0/0/0) 

(1,-1,0, 0, 0) 

(1,-1,0, 0, 0) 

(1, -1,0, 0, 0) 

MULTIPLY 




V 4 (0,0, 1,-1, 0) 

(0,0, 0, -1, 1) 

(-1,0, 0,0 1) 

(1,-1,0, 0,0) 

^ADD — 




0 

C (2) 
ov' ' 

- 1 

c (3) 

ov' ' 
2 

Cov<^) 
- 1 


Fig. 21. Cross-Correlation of Vectors u and V 

The cross-correlation operation is linear . This means that if u 
\v are vectors and is scalar, then 


C = aC 

au, V u, V 


and 


C , = C + C 

u , %• + W U , V u , w 


To prove Eq. (59), we use Eq. (57) to write 


m 


m 


C^u V + ■ "5^ u(t’)v(k + n = ttC^^^,(k) 

jfl i’ = l 




1 , - 1 , 0 , 0 ) 


Cuv<5) 

0 


, v, and 


(59) 


(60) 


Similarly, Eq. (60) is proved by writing 


m 

^=1 

m 

= y\(n [v(?+k) + w(?+k)1 

tTl -* 

m m 

= ^ ^ u(n v(Hk) + ^ ^ u(Ow(F+k) 

C=1 

= C (k) + C (k) 
uv uw 

As an example, the reader may verify that if u = (1. -1, 0, 0, 0), 

V = (0, 0, 1, -1. 0) and w = (0. 0. 0. 1. -l)then 

= 0) = 3C^ 3 X (0. -1. 2. -1. 0) 

and 

c ^ = c + c 

u, v+w u, V u, w 

= (0. -1, 2. -1. 0) + (0. 0. -1. 2. -1) 

= (0. -1. 1. 1, -1) 

Example 

As an example we consider the linear FM waveform discussed in the pre- 
vious section. Assume that a short linear FM return echo (of say 2M + 1 samples, 
where M is determined by the range pulsewidth T and the sample period At), 

^7 


from a point target is embedded in a much longer vector v (of say length N) 
consisting of all the terrain echo samples from a single pulse. We denote the 
2M + 1 nonzero linear FM samples ns 

s(-M), s(-M + 1) s(-l). s(0). s(l) s(M) 


Thus for k = 0, ±1 ±M. s{k) represents a sample from the return in 

Eq. (44), To see how we may conveniently express (44) in terms of k, the sam- 
ple index, let us translate the time axis in (44) according to eq. (51). so that 
time y is 0 in the middle of the echo return (44), Set 


B = A exp \ j 



and 


o = 



(61) 


(62) 


whore 


At = §4 (6 3) 

is given by (54), Then we have, for k = 0, ±1 ±M. 

•) 

8(k) = Be"'’’'' (64) 

(An easy way to sec the validity of (i'4) is to set time y = At in (44), correspond- 
ing to k = I in (t>4), Botli b in ((>4) and a in (n2) are commonly called the chirp 


■'.8 


rate. For example, in (44) or (10). the rate of change of radian frequency with 


respect to time is 


<t>" (t) = a 


(cf. eqs. (12), (13)). 

We wish to cross-correlate v with another vector u of similar form. More 
specifically, u contains the conjugate terms to Eq. (ol), with B replaced by 
unity. Thus somewhere in v» is a set of terms 

7(-M). ?(-M + 1) ilM - 1). "slM) 

whore 

7(k) = fork = 0. ±l ±M ('’5) 

All tho remaining terms ol* w 0, (The veetor u is also ol length N, ) 

Kxeept when tlie s terms in v overlap with the s terms in \i, the eross 
correlation veetor cletuenls are 0. When there is overlap, we may choose the 
index origin so that the index k runs iroiu -M to >M: 

\ U 

k = -Nt 

At first sight, it would appear tlial (<'(«) gives the correct lorimila tor 'U 

tor all . . However, rev all tliat stk il is nou/ero only tor 


(k+0 = 0. ±1 For k - M - i’ . s(k+0 = 0. Thus the Mghest summa- 

tion index in (66) is not M but M-i\ and the correct formula for the finite dura- 
tion discrete linear FM cross-correlation is 

M-f 2 

k=-M 

for the s and ~a overlap region 

- 2M s (k’+l) < 2M 


and 


C 


vu 


(t’+l) = 0 


(67b) 


otherwise. 

We show next that Fq. (67) may be simplified to 




li sin 

2 b 1 

(m + 4 ) 

1 f - bi’“ 


sin(b( ) 


for - 2 M 1 : (t+ 1 ) ^ 2 M 


( 09 ) 


The reader uiuonccrncd with details may skip to Fc|. ( 75 ). 
Ohserve first that Fit. (i' 7 ) simplifies to 


M-f 


E 


He- 


i 2 b k t 



(70) 
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Factoring out we have 


X V 

k=-M 


j2bki’ 


The right hand geometric series sums to the term in brackets below, which 


shows 


C (2+1) = Be 
vu 


^2 ^.j2bM2 ^j2b(M-2)2 ^ ^j2b2 


= Be 


jbt’^ |e-j2bM2 _ ^-^j(2bM2+2b2-2b2^) 


1 - eJ 


From both numerator and de nominator we factor out e* to obtain 


C (i’+l) = Be^' 
vu 


2 r -j2b(M+^)C j2b(M+j)C _j2b2^ 

jbi e - e X e ^ 

p-jb2 . Jb2 


Factoring out e ^ from the numerator yields 


, -j2b(M+|)2+jbC^ j2b(M+hk’-jbl’ 

^ , , . , , ib^*" .jbi e - e 

C„,(^+l) = BcJ <• ■> -ib,’ JbJ 

e - e 


j 2b(M+~)2-bi’^ -j ^2b(M+j)i-b.’^j 


ibi -jbi’ 
e* - e • 
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so that finally 


B sin 




[zb (m + 


sin (bg) 


(75) 


as promised. 

Note that Eq. (75) has a peak magnitude of |B • (2M + 1)| when ^ = 0. 
(Actually Eq. (75) is undefined at f = 0, but Eq. (67) may be easily evaluated 
when f = 0. ) Also the first null of appears when I satisfies 

2b (m g = ±iT 

or 


b?^ - 2b |m + ijt’ ± :t = 0 
or 


By completing the square we find 



(76) 


corresponds to the first positive null. Note that when |b| is large, f — 0. In 
fact, for fairly large |bl , we may write 
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where 


b 


TT 



« 1 


Using a Taylor expansion (as illustrated for example in Eqs. (85), (86)) we 
obtain 



Eqs. (76) and (77) show that at the first null, 


M + 4- 



J. 



TT 



(77) 


or 


libi ( M - t) 


UbM 


(78) 


I'or vahies ot t' greater than 


\Zb\ (M - T 

small compared to its magnitude ior values of f less tiian - j 

\lh\ IM " 


the magnitude of Kq. (75) is 

, Thais il 


we plot tlie luiu lion *l)j \crsus i , we obtain ttie response ol 1- ig . 

We see from I'ig. 11 that the e ross-eor relation of the return echo with 
the eoniugate esponential is the key to retries ine our lost range resolution. 

Tor. while' our transmitted ehirp was • 1 samples wide, tlie result ut eioss- 
eorrelating is onlv 1 ^ j p or j r samples wide 


essentially. 



ABOUT 

3dB 


2b (M + 1/2) 


2b (M + 1/2) 


Fig. 22, Plot of Square of Absolute Value of Linear FM 
Cross-Correlation Function 


In fact, the width of the half-amplitude (3 dB down) part of the curve in Fig. 22 
(shown with the dovible arrow - — ► ) is only [ 2 ' b ' \4 | ■ sppi'oximately. We define 
tills "half-width’’ loosely as resolution. More specifically, Fig. 22 shows the 
one-dimensional response to a point target. We can replace the abscissa units 
’’k”' in Fig. 22 by any convenient abscissa units desired, such as time t: 


t = TAt 


or slant range r: 


T 


If 


wo use tho t-axis, the half-width in t^tmits ie called the time- resolution. 
The half-width in r^-units is called the slant 
say that the dimension of the half-width 
When we recall that 


range spatial resolution, W'e also 
corresponds to one resolution cpH 


b = - f(At)^ 


Et,. ,62). then Et,. ,78, assutnes a special iornt tor anaiop ilnear FM 

signals. For, if *e let be the number of output samples within a range 
resolution cell, then 


rr 

2bM 




When tve realise that samples correspond to a time interval of f ^t 
we conclude that ^ 


seconds, 


TT-At 

2bM 


t 


r 


At 


IS an 


expresston for the time resolution of the radar. Sine. 


Ib| = j (At)^ 


we have 


(79) 


Time resolution 




Time resolution ' 

2M(At)a 


(80) 
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Now ZM(At) is the approximate linear FM pulse duration ^nd a is the radian 
frequency sweep rate, so that (2M(At)a) is the pulse radian bandwidth. Thus 

ig chirp bandwidth in Hz, and the time resolution (80) is the 

reciprocal 1/BW of the chirp bandwidth in Hz! Now we can see why the term 
’’compression” is used to describe the radar cross-correlation processing. 
After cross-correlation, the signal looks like Fig. 22 and has essentially a 
mainlobe time width of j= 1/BW. However, the original pulse time duration was 
T. Dividing the original pulse duration by the "processed” pulse duration gives 
the "compression ratio” CR: 

CR = t/ (ly/ Bw) = T X BW (81) 

Thus the compression ratio is the time -bandwidth product. 

We summarize: We send out a linear FM chirp of Hertz bandwidth BW. 
We wait for the echo which we then correlate with a conjugate replica of the 
original chirp, embedded in a vector of length equal to that of the echo vector. 

If a point target is in the reflecting field, it will produce a lobe structure as in 
Fig, 22, witli time resolution . The height of the mainlobe in Fig. 22 will 
be proportional to the reflectivity of the point target (amplitude factor B). The 
actual peak placement along the abscissa (labelled "0" in Fig, 22) depends 
directly on the slant range R to the target at the time of pulse transmission. 

If several point targets are found at varying locations along the image swath, 
linearity of the cross -correlation fvinction shows they will produce several 
separate mainlobes, appearing at several different places in the output cross- 
correlation vector, as in Fig. 23. 
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Fig. 23, Cross-Correlation Output When Two Point Targets Are 
Found at Different Slant Ranges 


Note froni Fqs. (ol) and (nP) that prior to taking absolute values, the 
cur\es in Fig. 2 3 all have azinuitli phase factors of the form 


exp 




I 

1 


exp 


1 4-r ) 

r ' I 


where R speeifies the target slant range (see Lq. (75)). And, in tact, absolute 
values are not taken right after range correlation. If tliey were, the all- 
important azitiuith pliase factor would Ije lost and the azimuth resolution would 
l)e limited by the antenna azimuth lieamwidtli. as in I'ig, 1 , However, the only 
reason w e took absolute values in the above explanation was to gi\ e the reader 
a picture (I'ig, 1 ') of what the range processed signal looks like. 



Wo ai*o almost ready to l)o)^in tlu' disoussion ol tbe a/dniuth proocssing. 
IIo\\e\cr, it will ho useful to lirst e haiMo terize the range' correlated return 
echo a bit more fully. Suppose' there is only the single point target PT in tiie 
beam wdicn our single pulse is transmitted. Then if f is the actual index of the 
cross-correlated return eclio, the actual form of tlie cross-correlation result 
is not Eq. (69) but is rather a translated version of Eq. (69), in wdiich t is 
replac ed by (t - q). We may w^rite 


C v^^( < ) = sin 


[zb (m t T)(t' - q) - b(,- - q)^ 
sin i b f i' - q] j 


(82) 


The translation index q (related to p in (51)) reflects the actual slant 
ranye to I'T, and is thus a fvinction of R, Recall also that 


B = A exp 


1 . 4t:R 

r' ^ I 


is also a lunction oi R, Ihvis tlu* location of the peak in I'ig. 24 is a function of 
R. and so is the azimuth phase tliat is attaelied to every point on the 

curve in Fig. 24. Tlte azimutli pliase is independent of iiulex i (ef. Fc|. (82)). 

1' rom the last remark it may be evident that optimum azimuth processing 
should make use ot every point on tlie ctirve ol Fig, 24, for e\'cry pulse echo 
in wliicli ]’l returns a signal. Sucli processing is currently regarded as too 
(. ompli'x, howe\er, and only the c|tli point is customarily used lor the azimuth 
processing. We will return to this point later (Fig, 27). 





IV, The Azimuth Processing 

The range processing involved operations on the information within.a single 

pulse echo. The azimuth processing will combine the information from a number 

j / 2 vj^R\| 

of echoes. The azimuth phasefactor = exp | j — c — /j 

key to this part of the processing. Recall that this factor arose from the mix- 
ing operation on the return echo. Although the azimuth processing is difficult 
to implement on a computer, it is conceptually simple, especially after the 
range processing has been understood. 

The first thing we have to do is examine where the data from a point 
target is found within a sequence of range correlated echoes. For simplicity, 
we assume the radar is mounted on an aircraft traveling past a point target PT 
on the ground. (The more complicated case of nonlinear or orbital motion is 
easily understood if the aircraft scenario is understood.) In Fig. 25 the slant 


TIMES OF PULSE 
TRANSMISSION 



Fig. 25, Slant Ranges as a Function of Aircraft 
Position Relative to Target 
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I'angcs to tlic target arc indexed as a function of the various positions 
successively occupied by the moving aircraft. We assume that only during 
pulses (-M) to (+M) docs the beam illuminate PT, Thus there are (2M+1) pulses 
carrying information about PT, We have labeled one slant range specifically 
because it represents the point of closest approach. The important thing to 
note is that the slant range R to PT changes with each pulse. We can express 
this quantitatively as follows. 

In Fig. 26, choose the x-axis along the flight path so that x = 0 when the 
aircraft is at the point of closest approach to PT. Then for any A/C position x, 
we have by the Pythagorean theorem that the slant range R to PT is given by 


R = 




{»') 



1 ig. 2i’. Rcl.iliontihip oi A t I’otiition to Slant Range 


'1 


In a practical SAR, x is very small compared to because we only need 
to consider values of x when PT is illuminated by the narrow azimuth beam. 
Therefore 


R. 


and we may write Eq, (83) as 


R = Ro Jl + 


(ro) 


(84) 


Eq. (84) has a radical of the form 


F(y) = JT+ y 


(85) 


where | y| !• Thus we may expand F(y) = Jl + y in a Taylor series about 
y = 0. Taking only the first two terms we have 


■t(y) ' 1 + ^ 


( 86 ) 


as a good approximation when |yl 1. 

2 

Substituting y = have 


R - R 




or 


R ' R * 

I' 2R,, 


( 87 ) 


\\ liii h ci\ os slant ranco R as a f\nu tion of A C azimutli position \. 


In l'<|« (S7). till' Niiricililo x is liiii'Cirly rdciti'd to tlic pulso index ni. whoio 
-M • 111 • M, riu' pulse index ni is in turn assoeiaied with tlie slant ranee 
in I'ip. 25. Wlien the aireralt is at position x, tlie slant ran.ee .eiven hy I'q. (87) 
will determine the loi ation ol tlie ran.ee eorri lation peak in l-ig, 24. In partii - 
ular. i| in (82) will he determined hy R whieh is itsi'll a c|u.idratii t'uni tion oi* x. 
Thus ii w e examine the sequenee ot ran.ee eorrelated echos associated with PT, 
we see that tlie location oi the peak in Fi.e. 24 will vary with the pulse index, as 
in Fi.e. 27 (wliere M = M. I'i.eure 27 represents tlie display of ranee correlated 
data over a sequence of pulses durine which PT is illuminated. Fach echo in 
Fig. 27 is a copy of the sinusoidal waveform of I2q. (75). where the location of 
the peak in liq. (75) is now a function of x (i. e. . the pulse index). .-\s mentioned 
earlier, (e. e. p. 7 .ind p. 5(1). we wish to >. omi'ine tlie intiirmation trom all 
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wish to choose the strongest possible signals from PT in each pulse. To 
obtain the strongest signal from PT. tlius, we must take the dotted samples in 
Fig. 27. These dotted samples follow an approximate parabola, as shown by 
Eq. (87). It is these dotted samples which we will extract and use for the 
azimuth correlation of PT. Associated with each dotted sample is an azimuth 
phase factor of the form: 


B = A exp 


1 j 4 ttR ) 

r 1 


( 88 ) 


(Recall from Eq. (75) and the subsequent remarks that each dotted sample has 
also a magnitude factor of (2M + 1. ) 

R is. remember, a function of azimuth position x. Thus, placing the dots 
of Fig. 27 in an "azimuth" vector 


Az = (Az(1). Az(2) Az(N)) 

where 


N = 2M f 1 


we may again examine the effects of cros 
phase factor is the conjugate to Eq. (88), 


s-correlating Az witli a vector Az whose 
i. c, , at position x, a has phase factor 


exp f ) 


1 , 4rrR ) 
^ 1 


(c-f. Eq. (81)). 

If we recall that 


R ' R„ 
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from Eq, (87), then we see that the correlation function takes the form 


exp 




±L 

\ 



( 89 ) 


Let us throw the correlation function into a more "digital" appearing form. 
Since presumably aircraft position x is related to velocity v and time t by 


X = vt 


(90) 


(provided t = 0 corresponds to the point of closest approach), we have for 
Eq. (89) 


oxp 


. 4it 




(^n) 


Let 


PRF = pulse repetition frecivieiuy 


so that 

At = = pulse repetition period 


is tlie lengtli ol time between pulses. (Rec all tiuit in sei tion III. At denoted the 
rant;. e samplin,i; period; here howe\er we use it to denote the a/imvitli sampling 

period.) I'hen it k = ±1, ±2. ±' ±M. the a/.imutl> releretue luiu tion ot 

Ikl. t"l) is gi\en as a luiu tion ot discrete multiples ot time At l>y 


A/ (k) 


exp 



\ " (kAt) 
JR 

t 1 





In Eq, (92), the constant phase factor 

I • 9t: ,, 

exp j J — ^^0 

is us\ially omitted. 

Now we can see that the same correlation process as described in Sec- 
tion III may be used. In short, the vector in Eq. ( 92 ) is embedded in a long 
vector, padded with zeros, called the azimuth correlation function. The con- 
jugate (received) signal, Az(k), given by 


Az(k) = A exp 




v'(kAt)^ \| 

/) 


for k = 0, ±1 ±M 


(93) 


is embedded in an ecjually long zero-padded vector. Note tluit the Az(k) are the 
dots in £'ig. 27. 

1 

If wc ol)ser\ e that Az(k) also has a constaiU phc\se factor exp | - j — — | 
which may be removed, then we sec that the essential cross-correlation is 
between 


A/*.(k) - i‘xp 




4-* 

-- ^ 


'lAtT 


iR 


(94) 


and 


i ’ ’ M 

A/.(io ^ exp — X - -ri^- ^ 


• 2 » 

I'.c], ('•?) is ol the lorm cxp'| il>k j where 


( 95 ) 


4- 


(At)"" 


Thus uc may use T<|. (uU) to eomlude that the result of aaimutl, eorrelation 


has the form 


I( ,’) = D sin 


4itx- (At) 

\R.t 


(’'' 4 - 2 ) 




( 90 ) 


where now the factor B is just the amplitvide factor .A.. From the results of 
Section III, we know the time resolution of Eq. (*50) is the reciprocal of the 
azimuth bandwidth. Multiplying this time resolution by v. the aircraft velocity. 

^ives tHo iiziniutH spntinl resolution. 

More specifically, we know the azimuth doppler chirp rate in Eq. is 

the second derivative of phase. Since azimuth phase <t>^ is 


> ? 

*a ■ \R. 


the chirp rate is 


> 



If the radar illuminates PT over a time T,^. 


the radian bandwidth is 


RBW^^ 


4’tv'" 


X T 


a 


( 97 ) 


aiul the llert;* bandwidth is RBW ^sl- or 


'"'a=tir ' To 




The time T is called the coherent integration time. It is related to the physical 
a 

distance L over which the aircraft travels while illuminating PT, by the simple 
a 

fo rmula 


(The length L is called the synthetic aperture length. ) Therefore, Eq. (98) 
may be rewritten as 


B\V 


a 


2vL 
a 


(99) 


Assuniing the antenna beam pointing direction is fixed, the maximum possible 

length L. over which PT is illuminated is the same as the width of the azimuth 

a 

beam on the ground. This width, from Fig. 2, is given as Thus for the 

maximum value of 1. » Fq. (99) becomes 

a 


BW 


a 


='■ “o 


i.et L now denote the physical antenna azimut) dimension. Since 1^^.^ = y . 
(cf. Fiq. (2o)). we have 


BW = ^ . 
a I. 


( 100 ) 


Recalling that tinic resolution is given as . we have then 


At = 4 ^ 

a 2v 


( 101 ) 


as ttie azimuth time resolution, and vAt . or 



as the best obtainable azimuth spatial resolution. L is the physical antenna 
azimuth dimension. 

If. on the other hand, we use a shorter length for L. than the maximum 
possible synthetic aperture length, then Eq. (99) implies a time resolution 
value of 


"^^^a ■ 2vL 


and an azimuth spatial resolution of vAt^ or 


Ax 



(103) 


(104) 


as a function of synthetic aperture length L^. The reader will note from 
Eq. ^104) that the best azimuth resolution is obtained with the longest synthetic 
aperture length L^. When the longest possible synthetic aperture is used. 

Eq. (102) applies and the resolution is liniited by the physical antenna azimuth 
dimension L. Therefore, in a S.-\R the smallest physical antenna leads to the 
best resolution. .-\s an example, the SEAS.AT spaceborne S.AR physical antenna 
is 10 meters along the azimuth. The synthetic aperture length, however, is 
more than 4500 meters. 

,\ few remarks .ire necessary licrc to clarify tlie simplifying assumptions 
leading to the azimutli impulse response gis en i>y l-'q. ('■'n), l.c|. ('h.) was 

deriv ed assuming a point target IM' is alone in the image tield and produces a 
linear IM signal emhedded in a longer vector ol zeros, .‘^ince eav li point target 
.done itic a/imutli piauim es a distiiut nvnilinear traiiv (like tlic lisliliones oi 
1 ig. Jh) it is evident tli.it tl.e .u lual i ro s s -c o r r elation process must l>e per- 
formcil with ».ire; for ex.imjvle. tU'.iders lamiliar witli II 1 v r os s -v or relation 


SLANT RANGE R 



P'ig. 28. Fishbone Data Traces of Successive .A,zimuth Point 
Targets at the Same Slant Range 


I 
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'J 

Q 

e( 
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techtiiques will recognize the impossibility of performing a fast correlation via 
an I'FT for the azimuth processing because the ’'fishbone*' data traces in 
Fig. 28 arc not colinoar. If the data samples of each "fishbone" arc used 
separately to perform a time domain correlation, however, the absence of 
interfering data from the neighboring fishbones will actually result in a "better" 
impulse response tlian that depicted in Eq. (9^). 

We finally observe that the two -ditnensional point target res]^onse 
r'itupulse resj’K'nse" ) of a synthetic a])erture radar is tlie product of the 
one-dimensional res])onses oi iT'M and Tlie index f is a range index 

in (T'O: in f is an a.‘imutli index. It is a useful ext'rcisi* for llie reader 

to try to \ isuali/i' tins iiiqndse vc:iponsr , 


f .0 




u 


T 


V. Some Additional Information 

In all the above analysis -se assumed a single point target PI at slant 
range R from the radar. The linearity o. cross-correlation in both the range 
and arlmuth processing ensures that at every position in the output array, the 
pixel value obtained is "reasonably" proportional to the reflectivity o£ the 
corresponding location on the terrain. The illustration of this effect txas 
already given in Fig. 23, where the two mainlobes in the figure (together with 
the sidelobes) add together to produce a curve which appears about the same 

as that in the figure. 

Additional simplifying assumptions we have used may become apparent to 
the reader as he continues to explore the field of SAR more deeply. However, 
it is hoped that arnred with the basic understanding provided by the above analy- 

Sis the reader will easily grasp the attendant subtleties. 

For exanrple. although an orbiting apacccraf, .SAR docs not permit the 
use of the simple geometry of Fig. 2b. the azimuth correlation tunctioi 
theless still takes a form si.nilar to (,1,. because the quadratic approximation 
is still useful. Thus the genera, SAR arimu.h correlation (unction has the form 

1 


Az(t) = exp 


' I ■ 


ith chirp rate in Hertz/ 


where f is the doppler cent-r frequency and f is the azimut 
second, some interesting probie.ns related to SAR processing are deserihed in 

the next section. 
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1 





VI. 


Some SAR ProcessinR Problems 


A, Speckle 

Because a SAR is coherent (i, e. , preserves phase information), the 
typical SAR image is marred with irregular variations in brightness which look 
like a laser speckle pattern on the wall. The cause is the same for a SAR and 
a laser. 

Each resolution cell in a SAR image (i. e. , each "pixel" or correlated 
output point) may be thought of as corresponding to a patch on the ground of 
finite area. Within this patch there are many "point" targets, each contributing 
to the overall reflectivity of the scene, as in Fig. 29. Suppose there are N of 



Fig. 29. C miposition of a Resolution Coll 

t>2 


PT PT PTm where N->1. For each k. 1 i k < N. returns 

them, PT p ^ N 

a signal to the radar of the form 

Si^(t) = cos^«})( t - + 6u) 


where 


R = nominal slant range to resolution cell 
*(,) = th,= normal oxpected phaaa function, aa in Section II. For 

example 


4>(t) 




0 ^ = Phaae adjuatmcnt. a real nuntbor related to tbe actual slant 

range ol PT^ owing to surface roughness, viewing angle, etc. 
(Th..- phase .adjustment 0^ also depends on dielectric and 

rcflection-phase-inversion characteristics of PT^). Thus 

0 is random variable, uniform on [O. In] . an.l independen. 
k 

of 0, for k ^ 

A = |A ! = Iho amplitude or reflection efficiency of PTj^. 

k k’ 


After mixing and correlation, the return from PT, i» estimated as 


■■k ° ''k '■■’‘P { ' °kl 


(lOo) 


Uy the linearity of cross correlation, th 
resolution t. ell is 


c ovcr.ill correlated return from tlie 


N N 

^ H ’•k - Xa ''l' ^ ' 

k=l k=l 


( 107 ) 



W t' in.iy I'xpi’i'Stf I'.u'li i*j^ in Cartesian form l)y 


= •'^1, f.it' 




I . 


1. + i Q, 


( 108 ) 


whore Ij and Qj^ are real, rims lai. ( 107 ) booomos 


N 



k=l 


h - 1 


( 109 ) 


or 


S = I + iQ 


( 110 ) 


where 


1 



( 111 ) 


O 



(Hi) 


From lk|. ( 110 ), w e may .eo hai k to Uie polar representation; thus 


S :: I I 


lO A I'xp 



(11 ') 


W luT I' 
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A - . 1 “ ! L'‘- 

V 


(INa) 


.1 n<l 


(I 


tan 


IC 
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- -t 


(IMh) 


Since by £qs. (Ill), (112), I and Q are sums of many random variables, wc may 
assume by the central limit theorem (e, g. , [4]) that I and Q are identically 
Gaussian. Thus 


is Rayleigh and 






( 115 ) 




is exponential, with mean equal to standard deviation. Thus the intensity 

of the resolution cell is a random variable depending on look angle, slant range, 

2 

etc. The mean, E(A ), would presumably provide a better estimate of the 
"true" reflectivity of the resolution cell. However, we only have one sample 
value of A^ from our correlation processing, and the signal -to -noise ratio of 
the exponentially distributed random variable A is 0 dB. To improve the situa- 
tion, we may try to obtain more samples for A . The interested reader is 
referred to the excellent paper . [s], for further details. 


B. p'oeusinu 

Both the range and azimuth correlation and signal functions were of the 
form (Sections II, III) 


s(t) = exp j j j 

where wc choose "a" to match the return signal value. In practice, especially 
for tlie azimuth processing, it is difficult to know the true signal value of "a" 
because it depends on orbital information (cf. i;c|. (9:^), for example). Mis- 
match betAveen tlu* correlation value "a" and the signal value "a" results in 



degradation of the impulse response Eq. (96). There exist various feedback 
schemes to iteratively refine the value of "a" based on perceived image degra- 
dation, The interested reader is referred to [6] for further details on these 
"automatic focusing" ("autofocusing") procedures. 

C. Interpolation 

As in Fig. 27, it is necessary to take data values along a curve for the 
azimuth correlation. Rarely does a sampled grid of SAR data contain exactly 
the right points. The result of using samples offset from the peaks in Fig. 27 
is amplitude modulation of the signal resulting in increased sldelobe levels in 
the impulse response (Eq. (96)) and consequent image quality degradation. The 
interested reader may refer to [7], [8] for more information bn interpolation 
as a possible solution. 

D. Roundoff Error 

Digital processing is of course performed in a processor with finite word 
lengths (and all the attendant dynamic range and accuracy limitations). Further 
details and problem solutions are found in [l] and [ 9 ]. 

E. Sidelobes and Weighting 

Because of the form of the impulse response in Eq, (96), image sidelobes 
may sometimes be mistaken for targets. More often they will limit the radar 
calibration accuracy. One solution is to "weight" the correlation functions. The 
sidelobes are thus reduced but the impulse response mainlobe is broadened. For 
further background, see [lO], Chapter 7, and [1 1] , [12], 
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